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TRANSITION,  MINIMUM  CRITICAL,  MINIMUM  TRANSITION,  AND 
ROUGHNESS  REYNOLDS  NUMBERS,  FOR  SEVEN  BLUNT  BODIES 
OF  REVOLUTION  IN  FLIGHT  BETWEEN  MACH 
NUMBERS  OF  1.72  AND  15.1 


by 

Neal  Tetervin 


ABSTRACT:  In  the  present  investigation  results  of  the  stabilit 

theory  for  laminar  flow  were  compared  with  tiansition  locations 
determined  from  heat-transf er  distributions  obtained  by  pre¬ 
vious  investigators  for  seven  blunt  bodies  of  revolution  in 
supersonic  flight.  The  comparison  shows  that  when  transition 
occurred  it  took  place  even  though  the  boundary  la:, er  ws-s 
calculated  to  be  very  stable  with  respect  to  small  disturbances 
for  the  entire  region  between  the  stagnation  point  and  the 
transition  location. 

In  every  case  transition  occurred  at  a  larger  boundary-layer 
Reynolds  number  than  the  estimated  minimum  Reynolds  number 
for  transition.  Consequently,  no  contradiction  of  the  as¬ 
sumption  that  there  is  a  minimum  transition  Reynolds  number 
and  no  disagreement  with  the  results  of  the  method  for  esti¬ 
mating  this  Reynolds  number  is  found. 

Five  of  the  seven  cases  considered  contained  useful  transition 
data.  A  first  examination  of  these  five  sets  of  data  seems  to 
indicate  a  connection  between  the  boundarx -layer  Reynolds 
number  at  transition  and  the  maximum  roughness  Reynolds  number 
ahead  of  transition.  A  further  examination,  however,  shows 
that  the  scatter  of  these  data  is  too  large  to  conclude 
statistically  from  only  five  sets  of  data  that  a  connection 
really  exists. 

The  boundary-layer  transition  Reynolds  number  was  found  to  be 
Influenced  much  more  strongly  by  the  maximum  roughness  Reynolds 
number  ahead  of  the  transition  point  than  by  the  local  wall 
temperature  ratio  at  the  transition  point. 
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Transition,  Minlmun  Critical,  Minimum  Transition,  and 
Roughness  Reynolds  Numbers,  for  Seven  Blunt  Bodies  of 
Revolution  in  Flight  Between  Mach  Nximbers  of  1.72  and  15.1 


This  report  presents  the  results  of  an  investigation  of 
transition  from  laminar  to  turbulent  boundary- layer  flow 
on  seven  blunt  bodies  of  revolution  in  flight  in  the  Mach 
number  range  between  1.72  and  15.1.  Mr.  J.  R.  Katz  pro¬ 
grammed  the  required  computations  for  the  IBM  704  electronic 
computer . 
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ILLUSTRATIONS 

29°  Hemisphere-Cone  In  Flight  at  a  Mach  Number 
of  3.14  and  a  Reynolds  Number,  Reoo>  8.94  x  10^ 

(a)  Boundary-Layer  Reynolds  Numbers 

(b)  Local  Mach  Number,  Local  Wall  Temperature  Ratio, 
and  Local  Roughness  Reynolds  Number 

50°  Hemisphere-Cone  In  Flight  at  a  Mach  Number  of 
4.7  and  a  Reynolds  Number,  Rego ,  of  8.03  x  10^ 

(a)  Boundary-Layer  Reynolds  Numbers 

(b)  Local  Mach  Number,  Local  Wall  Temperature  Ratio, 
and  Local  Roughness  Reynolds  Number 

Hemisphere-Cylinder  In  Flight  at  a  Mach  Number  of 
5.5  and  a  Reynolds  Number,  Rego ,  of  9.75  X  106 

(a)  Boundary-Layer  Reynolds  Numbers 

(b)  Local  Mach  Number,  Local  Wall  Temperature  Ratio, 
and  Local  Roughness  Reynolds  Number 

"l/lOtb-Power"  Nose  Shape  In  Flight  at  a  Mach  Number 
of  6.68  and  a  Reynolds  Number,  Rego»  of  8.15  x  10° 

(a)  Boundary-Layer  Reynolds  Numbers 

(b)  Local  Mach  Number,  Local  Wall  Temperature  Ratio, 
and  Local  Roughness  Reynolds  Number 

Flat-Face  Cone-Cylinder  In  Flight  at  a  Mach  Number 
of  14.5  and  a  Reynolds  Number,  Reoo ,  of  .821  x  10° 

(a)  Boundary-Layer  Reynolds  Numbers 

(b)  Local  Mach  Number,  Local  Wall  Temperature  Ratio, 
and  Local  Roughness  Reynolds  Number 

Spherlcal-Segment-Nose  Cylinder  In  Flight  at  a  Mach 
Number  of  15.1  and  a  Reynolds  Number,  Regg ,  of 
.855  X  10® 

(a)  Boundary-Layer  Reynolds  Numbers 

(b)  Local  Mach  Number,  Local  Wall  Temperature  Ratio, 
and  Local  Roughness  Reynolds  Number 

Elliptical-Nose  Cylinder  in  Flight  at  a  Mach 
Number  of  13.29  and  a  Reynolds  Number,  Rego »  of 
4.69  X  10® 
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(a)  Boundary-Layer  Reynolds  Numbers,  Sections 
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Ratio,  and  Local  Roughness  Reynolds  Number 
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8  Comparison  Between  Transition  Reynolds  Number, 
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10  Transition  Reynolds  Number,  Reo^T*  Wall 

Temperature  Ratio  at  Transition^  C^w/^e)^ 
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Body  and  Flight  Parameters 

Body  Pressure  Distribution,  Wall  Temperature 
Distribution,  and  Radius  Distribution  Data 

Surface  Roughness  Data 

Calculated  Boundary-Layer  Parameters 
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SYMBOLS 


speed  of  sound 
constant  -  equation  (17) 
constant  -  equation  (53) 
constant  -  equation  (17) 

constant  In  Sutherland's  formula  for  the  viscosity 

skln-frlctlon  coefficient,  Cf  *  tw/pw^e 

pressure  coefficient,  (pe/p©  “  Poo^Po) 

specific  heat  at  constant  pressure 

constant  In  equation  (54) 

enthalpy 

ratio,  i6*/d) 

mechanical  equivalent  of  heat,  778  ft-lbs  per  BTU 
height  of  roughness 


equivalent  sand  roughness  height 
E-l)/^(  Kg-l) 

dimensionless  wall  shear  parameter, 

(^AieX^w/^e)  (  S  y)^ 


reference  length 


2()(k-1) 


Mach  number 


dimensionless  correlation  number, 
-C^^/«'w)C‘^e/dx)  ( t^/te )  ^ ) 
momentum  parameter,  2  [n(H„+2) 
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Na 

P 

P '  ( t  ^ 
Pr 

q 

r 

1 

R 

Rcl 

Rek 

Rag 

*^®0,c 

Reo,B 

Rex 

ROw 

RCflo 

St 

Sw 

t 

u 

Uk 

u 

X 


Nusselt  number  (see  eq .  (39)) 
static  pressure 

-(l/Me)dMe/dx 

Prandti  number 

heat  transfer  by  conduction_to  surface  per  unit 
area  per  unit  time  (c  ^y)w 

recovery  factor 

radius  of  cross  section  of  body  of  revolution 

H/L 

reference  Reynolds  number,  uyiC/v^ 

roughness  Reynolds  number,  uk^/vk 

momentum  thickness  Reynolds  number,  Ue^/Vg 

minimum  critical  Reynolds  number  of  stability 
theory,  based  on  momentum  thickness 

minimum  transition  Reynolds  number,  based  on 
momentum  thickness 

Keynoids  number,  Uex/ve 

Reynolds  number,  Uqx/vw 

Reynolds  number,  UqqTj/Vqq 

Stanton  number  (see  eqs.  (43)  and  (46)) 

surface  temperature  ratio  parameter,  (T^/To  -  1) 

temperature 

velocity  parallel  to  surface 

lf2V 

u/uM 

distance  along  surface,  measured  from  origin 
of  boundary  layer 
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X  x/E 

y  distance  measured  perpendicular  from  surface 

a  exponent  of  Prandtl  number  in  Reynolds  analogy 

parameter  (see  eq.  (40)) 

B  non-dimensional  pressure  gradient  parameter 

ratio  of  specific  heat  at  constant  pressure  to 
specific  heat  at  constant  volume 

7  full  boundary- layer  thickness 

Z*  boundary-layer  displacement  thickness, 

^(  l-^/pe^e  )  dy 


boundary-layer  momentum  thickness, 
^pu/pe^e ( 1-u/ue) dy 


K 

X 

M 

V 

P 

T 

4 


thermal  conductivity 

dynamic  viscosity 
kinematic  viscosity 
mass  density 
shear  stress,  p  ^^u/  ^y^ 


Subscripts 

a  measured  along  axis  of  symmetry 

e  local  value  at  outer  edge  of  boundary  layer 

E  effective  value 

f  value  for  dp/dx  -  0 

k  at  distance  TE  from  surface 
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at  standard  sea  level  atmosphere 

maximum  value 

at  stagnation  point 

for  zero  heat  transfer 

at  transition  point 

associated  transformed  quantity 

value  at  surface 

free-stream  conditions,  ahead  of  bow  shock  wave 

all  quantities  with  a  "bar”  are  dimensional,  all 
others  are  non-dimensional 
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INTRODUCTION 

It  Is  well  known  that  the  heat-transfer  rate  for  a  turbu¬ 
lent  boundary  layer  is  larger  than  for  a  laminar  one  and  that 
the  difference  widens  with  Increase  in  boundary  layer  Reynolds 
number.  This  difference  In  heat-transfer  rate  is  often  Important 
enough  to  justify  the  effort  to  provide  an  extremely  smooth  sur¬ 
face  In  order  to  increase  the  probability  of  extensive  regions 
of  laminar  flow.  When  all  other  factors  are  fixed,  the  extent 
of  laminar  flow  usually  Increases  with  decrease  In  surface 
roughness  If  the  extent  Is  less  than  that  for  a  perfectly 
smooth  surface.  Even  for  a  perfectly  smooth  surface,  however, 
the  position  at  which  the  flow  changes  from  laminar  to  turbu¬ 
lent  cannot  at  present  be  calculated  theoretically.  Two 
positions  related  to  the  transition  position  can,  however,  be 
estimated.  One  is  the  position  ahead  of  which  transition 
cannot  occur  when  the  boundary  layer  is  €ixposed  to  very  small 
disturbances  of  a  particular  type.  This  position  can  be  cal¬ 
culated  by  use  of  the  results  of  the  stability  theory  (ref. 

(1)).  The  other  and  more  forward  position  is  the  one  ahead 
of  which  transition  cannot  occur  when  the  boundary  layer  is 
highly  disturbed.  A  rough  estimate  of  this  position  can  be 
made  by  the  method  of  reference  (2). 

Because  all  information  that  can  lead  to  a  more  accurate 
prediction  of  the  transition  point  is  valuable,  those  two 
theoretically  determined  positions  are  compared  with  transition 
positions  determined  from  heat-transfer  distributions  obtained 
by  previous  investigators  for  seven  bodies  of  revolution  in 
flight  at  supersonic  speeds.  The  effect  of  roughness  on 
transition  is  also  considered. 

ANALYSIS 

Derivation  of  Differential  Equation  for  » 

In  order  to  estimate  the  location  at  which  the  boundary 
layer  first  becomes  unstable,  and  the  location  ahead  of  which 
transition  is  supposedly  impossible,  it  is  first  necessary  to 
calculate  a  number  of  laminar  boundary-layer  parameters.  A 
convenient  method  is  Cohen  and  Reshotko's  (ref.  (3)).  The 
method  was  developed  for  a  gas  which  i.s  both  thermally  and 
calorically  perfect  and  which  has  a  Prandtl  number  of  unity 
and  a  viscosity  that  varies  directly  with  the  temperature. 
Although  Cp  and V  are  consequently  both  constant,  their  numeri¬ 
cal  values^are  not  specified  in  the  Cohen  and  Reshotko  method. 

The  form  of  the  boundary-layer  momentum  equation  used  in 
the  present  analysis  is  obtained  from  the  herein  corrected  form 
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Of  equation  (B5)  of  reference  (3). 
to  depend  on  x  and  Is  ~ 


The  correct  form  allows  X 


r\ 

Me  / 

f  .  J 

(1) 


The  value  of  IT  that  appears  in  the  exponent  K  is  an  effective 
value.  Its  determination  is  discussed  in  the  section  entitled, 
"Calculation  Procedure."  From  equation  (1)  there  is  obtained 
the  differential  equation 


All  quantities  In  equation  (2)  are  non-dimensional.  From  the 
definition  (see  eq.  (34)  of  ref.  (3)), 


_ _L^  al  Uc  ^ _ \__ 

dx  Ug  ax 


and  the  relations 


and 


(3) 


(4) 


it  follows  that 

I  /  to  \ _ L 

^  (  tj  Me,  dlX 


Moreover,  the  correlation  number  n,  which  Is 


defined  as, 

(eq.  (22),  ref. 


(3)) 
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can  be  written  as 


n  = 


due 

d  X 


(6) 


The  ratio  (vq/vw)  placed  in  a  more  convenient  form  by 

noting  that  for_a  theri^lly  perfect  gas 

Pur  ^  'ie 

Jo  po 

Because  the  static  pressure  in  a  boundary  layer  is  independent 
of  the  distance  normal  to  the  surface,  it  follows  that 

h 

The  ratio  C^q/v^)  can  then  be  written  as 

^  ^  |/  I  (7) 


From  the  assumption  of  Isentroplc  flow  at  the  outer  edge  of 
the  boundary  layer  it  follows  that 


Vc-/ 


(8) 


The  viscosity  ratio  (pw/Po^  assumed  to  be  given  by  the 
relation 


where 


_ 

A^o 


(9) 

(see  eq.  (4)  of  ref.  (3)) 


(see  eq.  (5)  of  ref.  (3)) 


The  viscosity  at  the  wall  is  thus  calculated  by  the  Sutherland 
relation.  By  use  of  equations  (3)  and  (9),  equation  (7)  becomes 


3 

CONFIDENTIAL 


CONFIDENTIAL 
NOLTR  62-25 


'Hr 


1 

io  ) 

Z. 

[l*^: 

'Mfj 

1 

(10) 


When  the  relatlOD  (10)  is  substituted  into  equation  (6)  and 
relation  (4)  is  used,  the  result  is 


n  = 

which  can  also  be  written  as 


Z 

3-z^e 

VE'I 

n  = 


when  equation  (3)  is  used. 


(11) 


By  use  of  equations  (5)  and  (11)  the  quantity  (^©/^e^ 
that  occurs  in  equation  (1)  can  now  be  written  as 


_  feiKTiiy 


n 


f'(¥) 


X 


\  1  N/V  i 

1+ 


(12) 


where  3-Ve 


Now  Introduce  a  quantity  ♦  related  to  the  non-dimensional 
momentum  thickness,  9  7/  Rol,  by  the  definition 


/A 


(13) 


Then 


n  _ 


4) 


4)fAe 


(14) 
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Equation  (2)  then  becones 


i. 

ckx 


4> 


-J 


=  Nl  + 


L^i  ■  1 

dx 


- 


X  \K-, 


K 


.(li) 


(15) 


After  some  manipulation  and  the  use  of  equation  (4) ,  equation 
(15)  can  be  written  as 


cl4>  _  n1 

cl  X 


+  4 


'/  \  » 


L  i  K 


_  _  '  ol  A 

dx  >  Jx 


(16) 


Equation  (16)  is  the  boundary-layer  momentum  equation  In 
a  form  that  contains  the  Cohen  and  Reshotko  parameter  N.  All 
the  quantities  In  equation  (16)  except  4  and  N  are  obtainable 
from  the  given  data.  Then,  if  4  and  N  are  known  at  one  value 
of  X  the  value  of  4  can  be  found  at  a  slightly  larger  value  of 
X  by  use  of  equation  (16).  This  value  of  4,  together  with  the 
given  pressure  and  wall  temperature  distribution  and  equation 
(11),  determines  n.  Because  N  depends  only  on  n  and  T^/To, 
the  integration  of  equation  (16)  can  then  proceed.  Once  n  and 
4  are  known  all  other  boundary- layer  quantities  can  be  calcu¬ 
lated. 

Derivation  of  Integral  for  Determination  of  4  Near  Stagnation 
Point 


The  stagnation  point,  where  and  R  are  zero,  is  a  singular 
point  of  equation  (16).  When  a  numerical  step-by-step  solution 
of  equation  (16)  is  begun  at  the  stagnation  point  the  result  is 
usually  a  variation  of  4  with  x  that  is  highly  oscillatory  and 
diverging,  at  least  for  small  values  of  x.  One  way  to  avoid 
this  difficulty  is  to  calculate  4  from  an  integral  instead  of 
from  the  differential  equation,  equation  (16).  This  integral 
can  be  developed  by  noting  that  the  wall  temperature  distri¬ 
bution  is  a  symmetrical  function  of  x  for  a  body  of  revolution 
at  zero  angle  of  attack;  all  the  bodies  of  the  present  investi¬ 
gation  were  supposedly  at,  or  close  to,  zero  angle  of  attack. 
Because  the  wall  temperature  is  a  symmetrical  function  of  x 
about  the  stagnation  point  it  follows  that 
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When  the  wall  tenperature  Is  Independent  of  x,  Cohen  and 
Reshotko  (ref.  (3))  suggest  the  approxlaatlon 

hi—  A"*"  •  (17) 

this  approximation  allows  equation  (16)  to  be  integrated  In 
closed  form  and  so  results  In  the  desired  Integral.  For  a 
constant  wall  tenperature,  the  term  dX/dx  In  equation  (16) 

Is  zero  and  so  drops  out.  In  the  present  case  the  wall  tempera¬ 
ture  Is  assumed  to  be  constant  from  x  -  0  to  a  large  enough 
value  of  X,  say  xj^,  to  allow  the  step-by-step  Integration  of 
equation  (16)  to  be  started  at  xi  without  the  difficulties 
usually  found  when  the  integration  Is  begun  at  x  ■  0,  the 
singular  point  of  equation  (16). 


Upon  use  of  the  approximation  given  by  equation  (17)  and 
the  relation 

t  71 

'r\--9rr-7:7:yr^  ) 


which  results  when  equations  (11)  and  (13)  are  used,  equation 
(16)  becomes 


Equation  (19)  Is  a  linear  first-order  differential  equation  for 
After  Integration  the  result  Is 


(20) 


for  the  condition,  -  0  at  x  -  0.  In  the  present  calculations 
the  use  of  equation  T20)  from  x  -  0  to  the  value  of  x  at  which 
Me  Is  equal  to  .05  was  found  to  be  satisfactory.  At  the  first 
few  values  of  x  very  near  zero,  the  values  of  ^  computed  by 
equation  (20)  usually  do  not  form  a  smooth  sequence  of  values, 
but  this  Is  not  important  because  each  value  of  ^  Is  Independent 
of  Its  values  at  smaller  x.  Usually,  the  variation  of  with 
X  becomes  sufficiently  smooth  before  x  reaches  the  value  at 
which  Me  Is  equal  to  .05. 
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The  values  of  A  and  B  in  equation  (17)  are  found  by  aaklng 
the  line  given  by  equation  (17)  tangent  to  the  curve  N(n,Sw) 
at  the  stagnation  point  (see  fig.  4  of  ref.  (3)).  The  value 
of  B  in  equation  (17)  is  then  found  from  the  relation 


which  follows  from  equation  (17).  Moreover,  at  the  stagnation 
point  of  an  axisymmetric  body 


(22) 


(see  page  14  of  ref.  (3)).  Therefore,  from  equation  (17)  it 
follows  that 


(23) 


Adaptation  of  Cohen  and  Reshotko  Method  for  Calculation  by 
Electronic  Computer 


In  previous  work  at  the  Naval  Ordnance  Laboratory  the  Cohen 
and  Reshotko  method  had  been  adapted  so  that  calculations  could 
be  made  by  the  IBM  704  computer.  The  adaptation  consisted  in 
expressing  the  Cohen  and  Reshotko  parameters  ),  (CfRev/Nu)pj,.j^, 

and  (d/d)||«o  iu  analytic  fora  by  the  use  of  "least- 
square**  polynomials  in  the  vtr tables  n  and  S^.  In  the  present 
work  these  polynomials  were  adjusted  to  eliminate  slight 
discontinuities  in  }  and  in  (CfRew/Nu)p2._i  at  n  -  0.  Moreover, 

in  order  to  handle  cases ^ with  values  of  n  less  than  -.7  or 
so,  the  polynomials  for  \  ,  (CfRew/Mu)px-.i  and  (6/0)||.o  for 
n  <C  -  .3  were  replaced  by  straight  lines  that  had  the  same  slope 
and  value  at  n  *  -.3  as  given  by  the  polynomials.  For  the 
quantity  the  modification  consisted  in  the  use  of  the 

value  of  Htj.  at  n  -  -.5  for  n  -  .5. 

Value  of  np 

In  order  to  begin  the  computation  the  value  of  Oq  must 
be  known;  this  value  is  used  in  equation  (23) ,  and  also  to  get 
the  initial  value  of  4  and  so  of  (0'^~R^)o  use  of  equations 
(18)  and  (13).  The  present  calculations  are  confined  to  bodies 
of  revolution  at  zero  angle  of  attack.  For  such  bodies,  the 
value  of  the  non-dimensional  pressure  gradient  parameter  B  is 
1/2  at  the  stagnation  point.  Consequently,  the  correlation 
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number  n,  which,  in  general,  depends  on  both  B  and  Is, 
for  a  stagnation  point,  a  function  only  of  Sw  An  analytic 
expression  for  this  function  was  obtained  by  fitting  a  least- 
squares  polynomial  of  the  second  degree  to  the  values  of  n 
for  B  -  1/2  for  all  the  values  given  In  table  2  of  refer¬ 
ence  (3) .  This  analytic  expression  was  Incorporated  Into  the 
IBM  704  program. 

Separation  Criterion 


Although  the  calculation  of  the  separation  point  was  not 
an  objective  of  the  present  investigation  a  criterion  to 
indicate  the  occurrence  of  a  calculated  separation  point  was 
incorporated  into  the  calculation  procedure.  The  criterion 
was  obtained  by  noting  that  because  the  flows  under  consideration 
do  not  begin  at  a  separation  point  the  value  of  the  shear  parameter 
1  decreases  as  the  pressure  gradient  parameter  n  increases  (see 
fig.  2  of  ref.  (3)).  That  is,  the  friction  coefficient  decreases 
with  increase  in  adverse  pressure  gradient.  If  n  increases 
sufficiently  there  is  eventually  reached  a  value  of  I  such 
that  a  further  decrease  cannot  occur  unless  n  is  decreased. 

For  a  decrease  in  adverse  pressure  gradient  to  cause  a  decrease 
in  friction  coefficient  is,  however,  physically  unrealistic  for 
a  flow  that  does  not  begin  at  a  separation  point.  Consequently, 
only  the  upper  branch  of  the  curves  of  figure  2  of  reference 
(3)  seen  physically  allowable  for  the  present  computations. 

Because  the  correct  criterion  for  separation,  namely,  1-0, 
is  thus  not  attainable,  it  was  assumed  that  a  fair  estimate 
of  the  separation  point  can  be  obtained  by  assuming  separation 
to  occur  at  the  smallest  allowable  value  of  1  .  On  each  curve 
of  \  against  n  for  constant  Sw,  the  smallest  value  of  I  occurs 
at  the  largest  value  of  n.  These  maximum  values  of  n  depend 
only  on  S^;  the  analytic  form  of  the  separation  criterion  was 
therefore  obtained  by  fitting  a  least-squares  polynomial  of 
the  fourth  degree  to  these  maximum  values  of  n.  In  the  present 
computations  separation  was  not  indicated  in  any  of  the  seven 
cases. 


CALCULATION  PROCEDURE 
Given  Data  and  Flight  Parameters 


The  calculation  procedure  requires  that  the  shape  of  the 
body  R(x),  the  pressure  distribution  Pe/Po(x),  and  the  wall 
temperature  ratio  Tw/Tq(x) ,  be  given  either  in  analytic  or 
tabular  form.  In  addition,  all  the  quantities  given  in  table  1, 
except  were  needed.  Tbs  temperature  Tqq,  and  the 

density  in  the  free  stream  ahead  of  the  body  were  either 

part  of  toe  experimental  data  or  else  were  obtained  from 
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tabulated  properties  of  the  standard  atmosphere  (ref._(4)). 
The  temperature  ratio  ,  and  the  pressure  ratio  Pn/Pooi 

were  read  from  the  charts  of  reference  (5)  when  the  velocity 
Uoo  either  greater  than  or  only  slightly  less  than  7000 
^t/sec,  the  lower  limit  of  the  charts.  For  smaller  values  of 
Uqq ,  the  tables  of  reference  (6)  were  used. 


Reoo :  The  free-stream  Reynolds  number,  Reoo> 


from  its  definition 


U.^L 


was  computed 


The  velocity  u^q  and  the  kinematic  viscosity  v^  were  either 
given  explicity  or  else  were  obtained  by  calcu^tion  from 
other  given  quantities. 


Rsl:  The  reference  Reynolds  number,  Re^,  was  computed 
from  the  definition 


where 


and 


QO 


(24) 


The  ratio  Pq/poo  calculated  from  Sutherland's  formula 

c 

(c=  I  “iS.fe'-E. 


and  the  previously  obtained  value  of  to/Too>  The  effect  of 
dissociation  on  the  ratio  |1o/moo  given  by  Sutherland's  formula 
can  be  estimated  by  use  of  table  6  or  figure  8  of  reference 
(7) .  In  most  cases  of  the  present  investigation  the  correction 
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to  Sutherland's  formula  was  less  than  ten  percent.  Because 
this  correction  is  small  and  because  a  more  Inexact  formula 
than  Sutherland's  is  used  in  the  Cohen  and  Reshotko  method, 
namely,  equation  (4)  of  reference  (3),  the  small  correction 
for  dissociation  was  not  used.  When  the  velocity  Uqo  was 
either  greater  than  or  only  slightly  less  than  7000  ft/sec, 
the  ratio  Poo/Po  read  from  the  app£oprlate  chart  of 
reference  (5;.  For  smaller  values  of  this  ratio  was 

calculated  from  the  relation  for  a  perfect  gas,  namely. 


^ E;  In  the  present  calculations  the  departure  of  air  from 
perfect-gas  behavior  at  high  temperatures  is  partially 
accounted  for  by  allowing  the  value  of  V  to  have  a  value 
other  than  1.4  In  the  relations  involving  this  ratio.  The 
value  of  used  to  obtain  the  relations  between  various  flow 
quantities  behind  the  nose  shock  is  called  g  and  is  assumed 
to  be  independent  of  x  and  y.  For  values  of  Mqo  loss  than 
about  four,  the  value  of  X  £  was  taken  as  1.4.  For  larger 
values  of  Mgo*  the  value  of  X g  was  initially  estimated  by 
use  of  the  charts  of  reference  (5)  by  first  locating  the 
stagnation  point  conditions  on  the  Mollier  chart  of  refer¬ 
ence  (5)  and  then  proceeding  along  a  line  of  constant  entropy. 
The  values  of  log  p/pi  along  this  line  of  constant  entropy 
were  plotted  against  the  corresponding  values  of  loe  p/pj 
The  result  was  approximately  a  straight  line  with  slope  y  £. 


Later,  a  simpler  method  for  calculating  Vg  was  used.  In 


this  method  the  value  of  p  g  was  chosen  to  give  the  correct 
value  of  the  non-dimensional  velocity  gradient,  (due/dx),  at 
the  stagnation  point.  This  value  of  5  g  was  obtained  by  noting 
that  from  the  equation  of  motion  for  the  flow  outside 
boundary  layer, 

^  tt 


the 


e< 


as 


d  X 


it  follows,  after  use  of  L' Hospital's  rule,  that 


(25) 
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Moreover,  from  equations  (3)  and  (8)  a  different  expression 
can  be  obtained  for  (due/dx)o.  This  expression  contains  V  ^ 
and  Is, 


y 

\  "[pa  J 

(26) 


The  value  of  uy,  the  maximum  attainable  velocity.  Is  equal 

to  '|f2Tio7j  where  Hq  is  the  stagnation  enthalpy.  The  value  of 
is  now  defined  to  be  that  which  makes  the  value  of 
(due/dx)^  calculated  from  equation  (26)  equal  to  the  value 
calculated  from  equation  (25).  Thus,  after  equating  equations 
(26)  and  (25) ,  the  result  for  X  g  is 


Experience  showed  that  V  v  could  be  calculated  more  quickly 
by  use  of  equation  (27)  than  by  the  procedure  involving  the 
use  of  the  Molller  chart  of  reference  (5).  Moreover,  the 
value  of  £  obtained  by  use  of  equation  (27)  was  close  to 
that  calculated  by  the  use  of  the  Molller  chart. 


value  of  the  stagnation  enthalpy,  Hq,  was  calcu¬ 
lated  by  use  of  the  relation 


(28) 


The  value  of  Cp  was  taken  as  7.725  BTU/slug/deg  Ranklne  and 
00  taken  equal  to  1.4 

dM^/dx:  In  order  to  integrate  the  momentum  equation, 

equation  (16) ,  the  quantity  dM^/dx  is  needed.  By  differentiating 
equation  (3)  there  is  obtained  the  relation 
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In  some  cases  the  data  were_given  as  the  pressure  coefficient 
ratio  Cp/Cpo(x)  rather  than  as  Pe/Po(x) .  In  these  cases  Pg/Po 
was  obtained  from  Cp/CpQ  by  the  relation  that  follows  from 
the  definition  of  Cp,  namely, 


(30) 


(dMe/dx)^;  At  the  stagnation  point  dll^/dx  is  needed  in 

order  to  calculate  (9*^  Note,  however,  that  it  Is  not 

needed  there  for  the  Integration  of  equation  (16)  because 
this  equation  is  not  used  until  a  larger  value  of  x.  Because 

-  /Pe  \ 

Pe/po  is  a  symmetric  function  of  x,  the  derivative  tl  — )  /dx  is 

^Po  ' 

zero  at  X  -  0.  Consequently,  relation  (29)  is  indeterminate 
there.  In  order  to  obtain  an  expression  for  (dMie/dx)^, 
L'Hospital's  rule  was  applied  to  equation  (29)  with  the  result 


(31) 


d^/dx^Cpg/po)  The  value  of  d^/dx^Cp^/p^)  ^  or  of 

d2/dx2(Cp/Cpo)  Of  in  equation  (31)  was  obtained 

either  from  an  analytic  expression  for  Cp/CpQ  or  from  tabular 
data  obtained  by  smoothing  the  values  of  (Pe/Po)  read  from  a 
graph.  When  smoothed  tabular  data  were  used,  the  value  of 

d2/dx2(Cp/CpQ)  Q  was  obtained  by  one  of  two  methods.  In  the 
first,  values  of  Cp/Cp^  were  read  at  equal  Intervals  in  x 
near  x  *  0,  a  difference  table  constructed,  and  the  second 
derivative  obtained  by  use  of  Newton's  forward  difference 
formula,  namely, 
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(32) 


This  formula  Is  obtainable  by  differentiating  equation  (5), 
page  192  of  reference  (8) . 


When  this  method  was  used  for  the  "l/lOth-Power"  Nose 
Shape  and  for  the  Elliptical-Nose  Cylinder  the  values  of  both 
(dyReL)o  Qo  form  a  "smooth"  extension  of  their 

values  for  larger  x.  It  was  found,  however,  that  the  values 
of  (^y^eL,)o  ^o  fair  smoothly  into  their  values  at 
larger  x  when  the  Cp/Cp^  distribution  for  small  x  was  repre¬ 
sented  by  the  parabola 


^  =  1—  (33) 


The  value  of  "a"  was  calculated  for  the  "l/lOth-Power"  Nose 
Shape  by  making  the  parabola  given  by  equation  (33)  give  the 
smoothed  value  of  (Cp/Cp^)  at  x  •  ,1.  For  the  Elliptical-Nose 
Cylinder  the  corresponding  value  of  x  was  ,05  instead  of  .1. 

Smoothing  of  Data 


By  experience  it  was  found  that  in  order  to  obtain  suf¬ 
ficiently  smooth  calculated  distributions  of  the  heat  transfer 
j  and  some  of  the  other  boundary- layer  parameters  with  x,  the 
Pe/Po  data  had  to  be  "smooth."  When  p^/poCx)  was  given  in  the 
form  of  a  mathematical  expression,  the  smoothness  requirement 
was  satisfied.  When  Pe/^o(^)  given  in  the  form  of  a 

graph,  pe/po  read  at  convenient  intervals  in  x  to  produce 

a  table.  These  tabular  values  were  then  smoothed  by  a  five¬ 
fold  application  of  formula  (1)  on  page  276  of  reference  (8). 
This  is  a  five-point,  third-degree,  least-squares  smoothing 
formula.  Although  it  was  not  certain  that  the  smoothing  pro¬ 
cedure  was  necessary  it  was  also  applied  to  the  radius  distri¬ 
bution  R(x) ,  and  the  wall  temperature  distribution  tw(x).  The 
smoothing  procedure  often  produced  a  slight  change  in  the 
original  data  distribution. 

Formulas  for  Boundary-Layer  Quantities 

Re^:  One  of  the  boundary-layer  quantities  of  interest  is 
the  boundary-layer  momentum  thickness  Reynolds  number  Re^, 
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defined  as 


By  introducing 
written  as 


the  reference  Reynolds  number  Re^^,  Reo  can  be 


or, 

and 


after  use  of 
(8)  as 


the  perfect-gas  law  and  equations  (3) , 


(4), 


(34) 


The  value  of  Re^  in  the  present  analysis  was  calculated  by 
using  the  Sutherland  viscosity  formula  for  the  ratio  Po^^e* 

Rek/h^ReL^^  :  Often  of  interest  is  the  roughness  Reynolds 


number,  defined  as  — 

■e  - 

A  convenient  expression  for  Rek  for  small  values  of  Oi/i)  can 
be  obtained  by  expanding  the  velocity  u  and  the  kinematic 
viscosity  V  In  a  series  In  y  and  keeping  only  the  first  power 
of  y.  To  the  first  order  In  y,  the  expression  for  Rek  becomes 


K 


(35) 


The  expression  (35)  Is  expected^to  be  sufficiently  accurate 
when  ¥  is  not  much  larger  than  6.  In  the  present  Investigation 
the  measured  values  of  ¥  on  all  the  bodies  were  always  less 
than  (S/7) . 


When  the  relation  _ 


(36) 
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which  is  equation  (21)  of  reference  (3) ,  is  used  together^ 
with  the  perfect-gas  law  and  the  approximation  that  *  Pwi 
equation  (35)  can  be  written  as  _ 

'  -L  \3 


K 


©A 


(37) 


When  equations  (3),  (4),  and  (8)  are  used,  equation  (37) 
becomes 


The  computations  made  by  the  IBM  704  machine  gave  the 
quantity  Ro^/k2ReL3/2,  From  this  quantity,  the  roughness 
Reynolds  number  Re^  was  calculated  for  the  desired  value  of 
k. 


(Nu/x) ;  The  heat  transfer  by  conduction  per  unit  area  of 
surface  per  unit  time  is  given  by 


The  Musselt  number  is  defined  as 

—  _ 

in  reference  (3).  In  the  present  analysis  this  definition  is 
generalized  to 


(39) 


because  of  the  high  temperatures  in  some  of  the  cases  analyzed. 
In  order  to  calculate  Mu,  the  relation 


(40) 
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which  Is  equation  (38)  of  reference  (3),  is  used, 
of  the  definition  of  Cf,  of  Ray,  and  the  relation 


and  equation  (36) , 


M-ur 

it  follows  that 


Upon  use 


(41) 


When  equations  (4)  and  (41)  are  used,  equation  (40)  becoses 


The  value  of  the  Prandtl  number  was  read  fros  figure  11 
of  reference  (7)  for  the  average  temperature  and  pressure 
on  the  surface  of  the  body.  The  value  of  a  was  taken  as  .4, 
the  value  suggested  in  reference  (3). 


Also  of  Interest  is  the  Stanton  number,  Stoo>  defined 


00 


(43) 


The  Stanton  number  can  be  expressed  in  terms  of  the  Nusselt 
number  by  eliminating  q  from  relations  (39)  and  (43) .  Upon 
use  of  the  definition  of  the  Prandtl  number,  the  result  is 


5t 


00 


nIo.  1 

*^r.r  ^ 


When  equation  (9) 


is  used,  equation  (44)  becomes 


Si 


to 


/\lu^  A  /  \ 


(44) 


(45) 


The  value  of  Pr^  used  in  equation  (45)  was  the  same  as  that 
in  equation  (42) . 
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The  Stanton  number  can  also  be  based  on  local  gas  properties 
at  the  outer  edge  of  the  boundary  layer.  In  this  case  the 
Stanton  number  is  denoted  by  Ste  and  is  defined  as 


(46) 


_  _  By  eliminating  q  between  equations  (40)  and  (46),  using 
PgUg  to  form  Re0,  and  then  using  equation  (34)  for  Reg,  the 
result,  after  using  equation  (9),  is 


(47) 


Here,  again,  the  value  of  Pr^  is  the  same  as  in  equation  (42). 


wr 


ille 


The  heat  transfer  was  calculated  by  use 
en  as  _  v 


of  equation  (43) 


(48) 


The  ratio  ISe/^o  calculated  from  the  relation  for  the 
conservation  of  energy  outside  the  boundary  layer,  namely, 


I 


2Te 

The  gas  is  assumed  to  be  thermally  and  calorlcally  perfect 
with  a  constant  effective  va.lue  of  V  called  V  £.  For  such 
a  gas  the  relation  between  He  and  the  speed  of  sound,  ae>  is 


(49) 


When  this 
result  is 


relation  is  used  in  the 
he  _ ! - 


relation  for  He/Ho  the 


(50) 
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The  ratio  Hr/He  is  calculated  trom  the  relation 

A  =  i-t-r-ii 

he  -2.  hg. 


Where  r  is  the  recovery  factor,  taken  equal  toYPr^, 
This  relation  becomes 


he. 


(51) 


when  equation  (49)  is  used.  Upon  the  use  of  equations  (50) 
and  (51) ,  equation  (48)  becomes 


In  order  to  partially  account  for  the  fact  that  Hw/Kq  is  not 
necessarily  equal  to  T^/Tq,  without  Increasing  the  length  of 
the  computation,  the  ratio  ‘E«/Fq  was  approximated  by  the 
expression 


(53) 


An  expression  for  b  was  constructed  by  imposing  two  conditions. 
The  first  was  that  at  the  stagnation  point 

Therefore  b  is  given  the  value  CHwo/^o)/C^wo'^o)  when  ^w  •  "^Wo* 
The  second  condition  was  that 


where 


to 


Therefore  b  is  given  the  value  unity  when 
conditions  then  are 


b  = 


“for  ■^v*^o 
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and  _  _ 

t>-  1  -for 


For  other  values  of  b  is  assused  to  vary  linearly  with 
The  result  for  b  then  is 


(54) 


where 


Because  the  wall  temperature  usually  fell  with  increase 
In  distance  from  ^he  stagnation  point,  more  accurate  calcu¬ 
lated  values  for  q  would  have  resulted  If  b  had  been  made 
equal  to  unity  at  -  Too  Instead  of  at  -  Tp*  This  change 
In  b  would  have  Increased  In  equation  (52)  and  so  would 

have  reduced  the  calculated  values  of  q. 


In  every  case  the^approprlate  calculated  quantity,  namely, 
either  Stoo>  St^,  or  q,  was  compared  with  the  distribution 
along  the  surface  obtained  from  the  experimental  temperature 
data  and  presented  In  references  (9)  to  (16).  For  these  seven 
cases  the  calculated  and  experimentally  obtained  quantities 
were  close  enough  together  to  Indicate  that  the  method  used 
to  calculate  the  characteristics  of  the  laminar  boundary  layer 
Is  sufficiently  accurate  for  the  purposes  of  the  present 
investigation. 

6*/B,  6/0:  The  boundary- layer  thickness  ratios  6*/0  and 
6/9  were  also  calculated.  The  expression  for  6*/0  is 


and  the  expression  for  6/9  Is 


(55) 


(56) 


Equations  (55)  and  (56)  are  equations  (40)  and  (41)  respectively 
of  reference  (3)  with  IT  replaced  by 
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T^/p^u^^:  The  local  friction  eoofflclont,  Is 

often  of  Interest.  The  relation 


ft®;* 

Is  used  with  relations  (36) ,  (9) »  (8) ,  (4) ,  and  (3)  to  obtain 


(57) 


6:  In  order  to  astlsate  the  value  of  Ree  below  which  the 
stability  theory  (ref.  (1))  predicts  that  all  snail  two-dl- 
nenslonal  disturbances  decay,  the  tables  of  reference  (17) 
were  used.  In  order  to  use  these  tables  It  Is  necessary  to 
know  the  value  of  the  non-dlnenslonal  pressure  gradient  paraneter 
B.  The  coBputatlons  by  use  of  the  Cohen  and  Reshotko  nethod 
(ref.  (3))  result,  however,  In  the  non-dlnenslonal  pressure 
gradient  paraneter  n  Instead  of  B.  The  paraneter  B  Is  a 
function  of  n  and  the  non-dlnenslonal  wall  tenperature  paraneter 
S^.  This  function  la  shown  In  figure  4  of  reference  (17).  In 
the  present  Investigation  8  was  expressed  as  a  function  of  n 
and  Sv  by  a  least-squares  polyncwlal  of  the  fourth  degree  In 
both  n  and  and  the  values  of  B  were  calculated  fron  this 
polynonlal . 


Ree,c*  nlnlnun  critical  Reynolds  nunber,  Ree^c 

the  Reynolds  nunber  below  which  all  snail  two-dlnennlonal 
disturbances  are  danped.  Once  the  value  of  B,  of  and 
of  3^  was  known  at  a  value  of  x,  the  corresponding  value 
of  Ree^c  obtained  fron  the  tables  of  reference  (17). 

Because  the  values  of  B,  and  usually  did  not  correspond 
exactly  to  an  entry  In  the  tables  of  reference  (17)  It  van 
necessary  to  Interpolate.  A  sufficiently  accurate  nethod  of 
Interpolation  was  first  to  change  the  tables  to  tables  of 
logxoR*e,c  instea*:  of  Re^^c*  A  linear  Interpolation  pro¬ 
cedure  In  B,  Mq,  and  Sy  was  then  used  to  find  the  value  of 
Ree^c  ^or  given  valuss  of  8,  ll«,  and  S^.  This  interpolation 
was  part  of  the  calculation  routine  and  was  nade  by  the  IBM 
704  electronic  conputer. 

Ro6,n*  I'ho  nlnlnun  transition  Reynolds  niusber  Ree^a  Is  the 
value  of  Rea  below  which  transition  supposedly  cannot  begin*  The 
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value  of  Re^  m  v&s  estimated  for  each  value  of  x  by  first  con¬ 
structing  a  iable  of  Reg  f  against  Mg  from  the  curve  for 
Re^  m  f  insulated ’surface  given  in  figure  2  of  reference 

(2) ’and  then  interpolating  to  find  the  local  value  of  Reg  ^  f 
for  the  local  value  of  Mg.  For  reasons  given  in  the  "Discussion" 
and  in  reference  (2)  only  the  curve  for  an  Insulated  surface 
was  used  even  though  the  surfaces  were  not  Insulated.  The 
value  of  Reg^m  was  then  obtained  by  the  relation 


where  the  ratio  (Reo  B)/(Re6,m,f)  was  obtained  from  the  curve 
marked  "equation  (67)"  of  figure  5  of  reference  (2).  This 
curve  was  first  converted  to  a  table  of  (Red,B)/(Reg  )  for 
equal  increments  in  B  and  an  interpolation  was  made  io  find 
the  value  for  a  required  value  of  B. 


Method  of  Integration  of  Differential  Equations 

In  order  to  calculate  the  value  of  the  integral  in  equation 
(20),  the  integral,  called  I,  was  calculated  from  the  equation 


Me  _ _ 


(58) 


This  differential  equation  was  Integrated  by  a  procedure  that 
was  already  coded  on  the  IBM  704  electronic  computer  (ref.  (18)). 
The  procedure  consisted  in  using  the  fourth  order  Runge-Kutta 
method  for  the  first  two  steps  in  x  and  then  using  the  fourth 
order  Adams-Houlton  predictor-corrector  method.  The  same 
procedure  was  used  to  integrate  equation  (16)  from  the  first 
point  at  which  it  replaced  equation  (20)  to  the  end  of  the 
range  of  x.  For  the  present  computations  the  IBM  704  elec¬ 
tronic  computer  used  eight  significant  figures  in  all  the 
calculations  and  the  results  were  printed  out  to  five  sig¬ 
nificant  figures.  The  step  size  in  x  was  .001  for  all  the 
calculations.  By  trial,  this  Interval  was  found  to  be 
sufficiently  small  so  that  a  doubling  of  the  interval  in  a 
few  test  cases  affected,  at  most,  only  the  fifth  significant 
figure  in  one  or  two  of  the  boundary-layer  parameters  at  a 
few  values  of  x. 
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BODY  DATA  AND  CALCULATED  RESULTS  IN  TABULAR  FORM 

Table  2:  Table  2  indicates  how  the  pressure  distribution, 
the  wall  temperature  distribution,  and  the  radius  distribution 
were  obtained  for  each  of  the  seven  bodies. 

Table  3:  In  table  3  are  given  the  "surface  roughness" 
data  for  the  seven  bodies.  Also  given  are  the  roughness  heights 
that  were  used  in  the  calculation  of  the  roughness  Reynolds 
number  Re^. 

Table  4:  The  quantities  listed  in  table  4  were  calculated 
by  the  methods  described  in  the  sections  entitled,  "Analysis," 
and  "Calculation  Procedure."  The  included  range  of  x  is  equal 
to,  or  slightly  greater  than,  the  range  for  which  the  flow  was 
judged  to  be  laminar.  From  the  listed  quantities  any  other 
quantity  for  which  a  formula  is  given  in  the  section  entitled, 
"Calculation  Procedure,"  can  be  computed.  To  get  Htr  snd 
btr/^tr  use  in  equations  (55)  and  (56)  it  is  necessary 
to  use  figures  7  and  8,  respectively,  of  reference  (3). 


SOME  FREE-FLIGHT  DATA  AND  CALCULATED  STABILITY 
AND  TRANSITION  RESULTS 

In  figures  1  to  7  is  shown  the  calculated  variation  with 
X  of  Re0,  Ree  Re^,  and  at  the  largest  Mach  number  at 
which  data  were  available  for  each  of  the  seven  bodies.  Also 
shown  is  the  calculated  minimum  critical  Reynolds  number. 
Regret  Aud  the  measured  wall  temperature  distribution. 

290  Hemisphere-Cone;  In  figure  la  is  shown  the  variation 
of  the  boundary-layer  Reynolds  number  Re©  and  the  minimum 
transition  Reynolds  number  Rea^m  with  the  non-dimensional 
distance  x  from  the  stagnation  point  for  a  hemisphere-cone 
with  an  Included  angle  of  290  (ref.  (9)).  Also  shown  is  the 
range  of  values  of  Re©,c*  Transition  occurred  somewhere 
between  x  -  .6632  and  x  "  .7850,  that  is,  between  380  and  45® 
from  the  stagnation  point,  at  a  value  of  Re©  between  642  and 
742.  In  this  region  of  the  body.  Re©  ^  was  more  than  100 
times  as  large  as  Re©.  Ahead  of  the  iransltion  region,  this 
ratio  was  still  larger.  Transition  therefore  occurred  even 
though  the  boundary  layer  was  estimated  to  be  very  stable 
with  respect  to  the  small  two-dimensional  disturbances  used 
in  the  stability  theory.  It  Is  apparent  from  the  figure  that 
Re©  exceeded  Re©,ii  beyond  about  3.2®  from  the  stagnation 
point. 
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Id  figure  lb  Is  shown  the  variation  of  the  local  Mach 
number,  Me,  the  local  ratio  of  the  surface  to  stagnation 
temperature,  ^w'^ot  local  roughness  Reynolds  number 

Rej^.  In  the  transition  region.  Me  li^  between  about  .8  and 
1.0.  The  wall  temperature  ratio,  T|y/to»  varied  from  about 
.565  at  the  stagnation  point  to  about  .515  at  x  .7850.  The 
maximum  value  of  the  roughness  Reynolds  number  was  .0157  near 
X  *  .60;  this  roughness  Reynolds  number  is  based  on  a  roughness 
height  of  three  mlcrolnches  (see  table  3) . 


The  data  given  in  table  4  for  the  lower  Mach  numbers  for 
this  body  also  show  that  transition  occurred  where  Ree^c 
very  much  larger  than  Re^.  It  is  noted  that  at  Macn  numbers 
of  2.32  and  2.47  there  was  a  region  of  the  body,  roughly 
between  x  -  1.5  and  x  -  1.9,  in  which  Rea^c  less  than 
Re^.  In  both  of  these  cases  transition  occurred  further  back, 
where  Re^^c  was  larger  than  Rea. 


50<^  Hemisphere-Cone:  In  figure  2a  is  shown  the  variation 
of  Re^ ,  Re^  ,  and  Rea  ^  hemisphere-cone  with  an 

included  an^Ie  of  50°  (ref.  (10)).  Transition  began  somewhere 
between  x  -  .3379  and  x  -  .5222,  that  is,  between  19.3^  and 
29.9^  from  the  stagnation  point.  In  this  region,  Rea  lay 
between  303  and  436.  The  critical  Reynolds  number,  Rea,c» 
was  more  than  100  times  as  large  as  Rea  ahead  of,  and  in  the 
first  part  of  the  transition  region.  Transition  is  supposedly 
Impossible  ahead  of  about  3.30  from  the  stagnation  point,  the 
region  in  which  Rea  was  less  than  Re^  At  the  lowest  Mach 
number,  namely  at  2.5,  the  data  in  ta6le  4  indicate  that 
this  point  moves  back  to  about  4.3<’. 


In  figure  2b  is  shown  the  variation  of  M^ ,  tw/To,  and  Rek • 
The  local  Mach  number,  14,  in  the  region  in  which  transition 
began  was  between  .40  and  .63.  The  local  wall  temperature 
ratio  fell  slightly  from  .46  at  the  stagnation  point  and  then 
rose  to  .57  near  the  30^  station  where  the  flow  was  definitely 
in  transition.  The  maximum  value  of  Rek  was  equal  to  about 
11.1  and  occurred  at  about  x  -  .36,  just  beyond  the  last 
thermocouple  at  which  the  flow  definitely  was  laminar.  The 
values  of  Rek  were  calculated  by  use  of  a  roughness  height 
of  70.7  mlcrolnches.  This  number  was  obtained  by  multiplying 
the  measured  rms  values  of  25  microinches  by '{T?  in  order  to 
obtain  the  peak  height;  the  factor , 'jTs'.  follows  from  the 
assumption  that  the  roughness  can  be  approximated  by  a  "sine- 
wave"  shape  (ref.  (19)).  The  value,  25  mlcrolnches,  was 
obtained  by  use  of  a  profilometer  before  the  surface  was 
oxidized  to  stabilize  the  emissivity.  The  oxidation  probably 
increased  the  roughness  height  and  so  the  actual  maximum 
roughness  height  might  easily  have  been  larger  than  70.7,  the 
value  used  to  calculate  Rek. 
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Just  as  for  a  flight  Mach  number  of  4.7,  transition  at 
the  lower  Mach  numbers  occurred  between  x  ■  .3379  and  x  -  .5222 
(see  table  4).  In  the  region  in  which  transition  began,  the 
critical  Reynolds  number,  Red,ct  is  more  than  60  times  greater 
than  Ree  for  all  five  Mach  numbers  for  which  data  are  given  in 
table  4. 

Hemlsphere-Cy  Under 

In  figure  3a  is  shown  the  variation  of  Ree,  Ree^n,  and 
Ree^c  with  X  for  a  hemisphere-cylinder  (ref.  (11)).  The 
results  from  this  flight  do  not  allow  transition  to  be  more 
precisely  defined  than  to  say  that  it  occurred  somewhere 
between  the  stagnation  point  and  the  11  1/2^  location  (x  - 
.2008).  Consequently,  transition  occurred  at  a  value  of  Reg 
less  than  220.  Transition  occurred  even  though  Re^  c  was 
more  than  1700  times  as  great  as  Rea.  The  data  in  iable  4 
indicate  that  for  all  five  Mach  numbers,  Re^  q  was  at  least 
800  times  as  large  as  Rea*  Transition  is  supposedly  impossible 
ahead  of  about  2.8^  from  the  stagnation  point.  Transition 
therefore  occurred  not  sore  than  about  9^  behind  the  most 
forward  possible  point.  At  the  lowest  Mach  number,  namely 
3.0,  the  data  in  table  4  indicate  that  the  most  forward 
possible  transition  point  moved  back  to  about  4.3^  from  the 
stagnation  point. 

In  figure  3b  is  shown  the  distribution  of  M«,  ty/To,  and 
Re^.  Transition  occurred  in  a  region  in  which  the  local  Mach 
number  was  less  than  .24.  The  wall  temperature  ratio  was  equal 
to  .24  at  the  stagnation  point.  The  roughness  of  the  body 
was  stated  to  be  about  three  to  five  microinches  near  the 
stagnation  point  and  about  five  mlcroinches  further  back. 
Scratches  of  the  order  of  15  mlcroinches  existed  behind  the 
stagnation  region.  These  numbers  were  obtained  by  use  of  an 
interferometer  microscope,  not  on  the  body  itself,  but  on  a 
sample  of  the  body  metal  polished  in  the  same  way  as  the 
body.  In  reference  to  the  early  transition  in  this  test, 
the  original  investigators  (ref.  (11))  suggest  that  perhaps 
irregularities  other  than  the  five  microinch  surface  roughness 
caused  the  early  transition. 

"l/lOth-Power"  Nose  Shape 

In  figure  4a  is  shown  the  calculated  data  for  a  body, 

(ref.  (12)),  whose  nose  shape  is  defined  by  the  equation 
given  in  table  2.  At  x  -  .416  the  flow  was  laminar  but  at 
X  -  .611  transition  was  already  well  along  toward  completion. 
Transition  therefore  began  at  a  value  of  Ree  between  277  and 
417.  At  the  lower  Mach  numbers  transition  occurred  further 
back.  In  fact,  at  the  lowest  Mach  number  transition  had  not 
occurred  at  x  1.920,  the  location  of  the  rearmost  thermocouple. 
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When  transition  did  occur,  Re0  q  was  more  than  250  times 
larger  than  Re^j  for  the  entire  length  of  laminar  flow  for 
all  Ifach  numbers.  The  indication  from  figure  4a  is  that  the 
value  of  Re^  was  less  than  Rea^n  for  x  ^  .1.  At  the  lowest 
Mach  number,  namely,  1.72,  the  corresponding  region  was  that 
for  X  ^  .15.  The  information  in  figure  4b  indicates  that 
transition  occurred  between  a  local  Mach  number  of  about  .3 
and  .3.  The  wall  temperature  ratio  varied  from  about  .26 
at  the  stagnation  point  to  about  .40  at  x  -  .611.  The  maximum 
roughness  Reynolds  number  in  the  region  in  which  the  flow 
definitely  was  laminar  was  equal  to  1.3. 

Flat-Face  Cone-Cylinder 


In  figure  5a  is  shown  the  variation  of  Reg,  Reg  ^  and 
Reg  Q  with  X  for  a  body  consisting  of  a  flat  face  followed 
by  a  cone  which  is  followed  by  a  cylinder  (ref.  (13)).  The 
last  thermocouple  on  the  cone,  at  x  -  4.493,  indicated  laminar 
flow.  A  calculation  of  the  heat- transfer  rate  at  x  -  5.290 
under  the  assumption  that  the  flow  was  turbulent  results  in 
a  heat-transfer  rate  of  about  44  BTU/sq  ft/sec.  The  calculated 
value  for  laminar  flow  is  about  five  BTU/sq  ft/sec  and  the 
value  given  in  reference  (13)  is  about  33  BTU/sq  ft/sec.  Conse¬ 
quently  the  flow  not  only  was  not  laminar,  but  probably  was 
completel>  turbulent.  It  is  remarked  that  the  value  of  x  for 
station  (13)  in  reference  (13)  should  be  5.290  instead  of  4.77, 
the  value  given  in  reference  (13).  At  x  -  4.906  the  calcu¬ 
lated  value  is  about  44  BTU/sq  ft/sec  for  turbulent  flow, 
about  five  for  laminar  flow,  and  the  value  given  in  reference 
(13)  is  about  73.  Therefore,  although  the  flow  may  not  have 
been  completely  turbulent  at  x  -  4.906,  it  does  not  appear 
to  have  been  laminar.  Transition  therefore  occurred  between 
X  -  4.493  and  x  -  4.906.  Similar  calculations  for  flight  Mach 
numbers  of  13  and  10  lead  to  the  same  conclusion,  namely, 
that  transition  began  between  x  -  4.493  and  x  -  4.906. 

Transition  thus  occurred  near  the  intersection  of  the  cone 
and  the  cylinder,  the  region  in  which  Reg  increased  rapidly 
by  a  factor  that  varied  between  about  3.2  at  a  flight  Mach 
number  of  14.5  to  about  2.2  at  a  Mach  number  of  10.  Conse¬ 
quently  transition  occurred  at  a  value  of  Reg  that  was  in  the 
neighborhood  of  500  even  though  Reg^c  infinite.  Beyond 

the  first  43  percent  of  the  front  face,  Reg  exceeded  Reg^n- 

In  figure  5b  is  shown  the  distribution  of  M^ ,  Tw/Tq,  and 
Re)(.  The  maximum  value  of  Rejc  is  about  1.7  and  occurs  on 
the  radius  of  the  nose  corner.  A  smaller  local  maximum  occurs 
at  the  cone-cylinder  junction;  its  value  is  about  .23.  The 
wall  temperature  ratio  was  about  .13  near  the  stagnation 
point  and  fell  to  about  .10  at  the  end  of  the  conical  portion 
of  the  body . 
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Spherical-Segment  Nose 


In  figure  6a  is  shown  the  variation  of  Red,  Re^  and 
Re^  Q  with  X  for  a  body  made  up  of  a  spherical  segment  followed 
by  k  cylinder  (ref.  (14)).  The  flow  appeared  to  be  laminar 
over  the  spherical  face  and  over  the  cylinder  at  least  as  far 
back  as  X  -  .716,  the  location  of  the  rearmost  thermocouple. 

Note  that  Re^  was  small  over  the  entire  body.  The  maximum 
calculated  value  is  about  170  Just  beyond  the  face-cylinder 
Junction.  The  value  of  Reg  then  dropped  rapidly  to  a  minimum 
near  20  and  then  increased  slowly  to  the  rear.  The  values  of 
Re^  were  so  small  that  although  they  exceeded  Rea  a  at  about 
80  from  the  stagnation  point  they  were  never  more' than  about 
2  1/2  times  as  large  as  Rea  b*  Moreover,  over  the  cylindrical 
portion  of  the  body  Rea,m  exceeded  Re^.  It  is  noted  that  these 
values  of  Rea  q  are  calculated  for  an  Insulated  surface  and  so 
may  be  too  large  because  the  ratio  Tw/Tq  is  very  low  in  the 
present  case.  Even  if  the  values  of  Red,m  were  only  half  wf 
those  shown,  Rea  would  still  be  less  than  Rea  n  over  most  of 
the  cylindrical  portion  of  the  body  and  not  much  larger  over 
the  face.  In  the  present  case  the  boundary  layer  remained  laminar 
where  stability  theory  predicted  it  should. 

In  figure  6b  is  shown  the_dlstrlbution  of  M^,  a-^J 

Reic.  Note  the  low  values  of  this  ratio  was  not  greater 

than  .12  over  the  entire  body.  Some  data  indicate  that  when 
the  ratio  Tw/Tq  is  small,  of  the  order  of  .25  or  so,  transition 
occurs  far  forward  on  the  body  (refs.  (20)  and  (21)).  In  this 
and  in  the  previous  case,  where  the  wall  temperature  ratio  was 
also  near  .12,  laminar  flow  existed  in  spite  of  the  low  value 
of  this  ratio.  It  is  noted,  however,  that  although  a  large 
portion  of  the  body  was  covered  by  laminar  flow,  the  maximum 
values  of  Rea  for  the  laminar  flow  were  small.  Also  to  be 
noted  is  that  in  the  present  and  in  the  previous  case  the 
Reynolds  number  per  foot,  which,  in  reference  (21),  is  stated 
to  be  important  in  the  "transition-reversal"  problem,  was 
low. 

Elliptical-Nose  Cylinder 

In  figure  7a  is  shown  the  variation  of  Rea,  and  Re^  b 
with  X  for  a  body  with  an  elliptical  nose  followed  by  a’ 
cylinder.  This  body  is  called  "3-204-1"  in  references  (15) 
and  (16).  The  smallest  value  of  Reg  ^  over  the  region  of 
laminar  flow  is  also  listed.  For  this  body,  calculations 
were  made  for  two  meridian  sections.  Along  section  ABCEA* 
the  measured  surface  roughness  height  was  1/2  microinch, 
root-mean-square.  Along  section  ALM  the  measured  roughness 
height  was  1/2  microinch,  rms  up  to  x  ■■  .3491;  for  x  >  .3491 
the  roughness  was  45  microinches,  rms. 
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Along  section  ABCEA',  the  smooth  section,  transition  began 
near  x  •  1.815.  The  nose-cylinder  Junction  Is  at  x  ■■  1.115; 
consequently,  laminar  flow  covered  the  face  and  a  portion  of 
the  cylindrical  afterbody  equal  In  length  to  .7  of  the  cylinder 
radius.  Transition  began  at  a  value  of  Re6  of  about  1080. 

Here  again,  transition  began  In  a  region  In  which  Red,c  was 
very  much  larger  than  Re^.  Transition  was  supposedly  impossible 
ahead  of  x  .086,  the  region  In  which  Re^  <  Re^^ni. 

In  figure  7b  la  shown  the  distribution  of  local  Mach 
number  Me,  the  ratio  of  the  local  wall  temperature  to  the 
stagnation  temperature  local  roughness  Reynolds 

number,  Reic,  for  the  smooth  section,  section  ABCEA*.  The 
peak  In  the  distribution  of  M«,  at  x  -  1.2,  was  Introduced  by 
the  procedure  for  smoothing  the  pressure  distribution  data. 

A  small  change  In  the  pressure  distribution  near  the  face- 
cylinder  Junction  Is  not  believed  to  affect  significantly 
any  conclusions  which  depend  on  the  values  of  Reg  or  Rej^ 
further  downstream. 

The  values  of  were  low;  they  varied  from  about  .2 

near  the  stagnation  point  to  about  .13  near  x  -  1.858.  The 
largest  value  of  Re^  was  about  .0182  and  occurred  at  x  -  .955. 

In  the  region  In  which  transition  began  the  value  of  Re^  was 
about  .0004. 

In  figure  7c  Is  shown  the  distribution  of  Reii;  for  section 
ALM.  At  the  edge  of  the  roughness  patch,  at  x  .3491,  the 
roughness  Reynolds  number  Jumped  from  about  .0030  to  about 
24.5.  Transition  began  somewhere  between  x  -•  .4190  and 
X  ■■  .6290,  not  far  from  the  start  of  the  roughness  patch. 


DISCUSSION 

Because  the  boundary-layer  Reynolds  number  Is  zero  at  the 
stagnation  point,  the  flow  Is  always  laminar  there.  Moreover, 
any  disturbance  of  the  laminar  boundary  layer  will,  because 
of  the  smallness  of  the  Reynolds  number  near  the  stagnation 
point,  die  out  as  It  proceeds  downstream  with  the  flow. 
Consequently,  transition  cannot  occur  very  close  to  the 
stagnation  point.  Further  back,  however,  the  local  boundary- 
layer  Reynolds  number  usually  becomes  large  enough  for  dis¬ 
turbances  to  the  laminar  flow  to  either  cause  transition 
almost  Immediately  If  they  are  large  enough,  or  to  grow  as 
they  proceed  downstream.  The  disturbances  that  grow  as  they 
move  downstream  usually  result  In  transition  somewhere  down¬ 
stream  of  the  location  of  Initial  amplification. 
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Conparison  with  Results  of  Stability  Theory 

Although  the  location  where  very  small  disturbances  begin 
to  grow  can  be  calculated  by  the  stability  theory  (ref.  (1)), 
the  transition  position  cannot  be  calculated  because  the 
disturbances  soon  become  too  large  for  the  stability  theory 
to  apply.  The  stability  theory  is  at  present  Halted  to  wave- 
like  disturbances  whose  amplitude  is  small  enough  to  allow 
the  equations  of  motion  and  energy  to  be  linearized.  In 
the  present  Investigation  the  results  of  the  stability  theory 
are  used  to  estimate  the  distribution  of  the  critical  Reynolds 
number,  Ree^c*  along  the  body  and  thus  also  the  point  at  which 
Ree  exceeds’ Red, c*  critical  Reynolds  number  of  the  present 

analysis  is  the  Reynolds  number  below  which  all  small  two- 
dimensional  wave-like  disturbances  die  out.  Above  this 
Reynolds  number,  disturbances  with  the  proper  frequency  can 
grow.  Consequently,  if  Re^  exceeds  Ree  q,  disturbances  can 
grow  and  the  flow  is  unstable.  Althougfi  Dunn  and  Lin  (ref. 

(1))  have.  In  their  stability  theory  for  compressible  flow, 
also  treated  disturbances  that  travel  at  an  angle  to  the  main 
stream,  the  present  Investigation  estimates  Re^  q  only  for 
disturbances  that  travel  in  the  direction  of  the’main  flow. 

It  is  not  clear  Just  what  sort  of  a  spiral  path  a  disturbance 
traveling  at  an  angle  to  the  main  flow  would  take  on  a  body 
of  revolution. 

Although  the  stability  theory  Is  derived  for  the  flow  over 
an  Infinite  plane  It  is  applied  in  the  present  investigation 
to  the  flow  over  a  body  of  revolution.  About  20  years  ago, 
Pretsch  showed  that  if  the  boundary-layer  thickness  over  a 
body  of  revolution  is  a  small  fraction  of  the  radius  of 
curvature  of  the  surface,  the  stability  theory  gives  the 
same  critical  Reynolds  number  for  the  flow  over  a  body  of 
revolution  as  for  the  flow  over  an  infinite  plane  (ref.  (22)). 
Pretsch  derived  his  result  for  incompressible  flow.  It  is 
assumed  that  the  same  result  is  valid  for  compressible  flow. 

In  the  present  analysis  the  values  of  Re^  ^  were  obtained 
from  the  tables  of  reference  (17).  These  values  were  calcu¬ 
lated  by  use  of  Lees'  formula  for  the  rapid  estimation  of 
Re0,c  (!))•  Lees'  formula  is  approximate  and  is  probably 

a  useful  approximation  only  when  the  local  Mach  number  outside 
the  boundary  layer  is  less  than  or  not  much  larger  than  unity. 
It  is  remarked  that  for  Mach  numbers  greater  than  about  two, 
even  the  exact  method  of  calculation  of  Re^  ^  appears  still 
to  be  under  development.  For  certain  condlilons  Lees'  formula 
predicts  that  Rea  q  Is  infinite.  It  turns  out  that  the  formula 
provides  a  good  estimate  of  these  conditions  (ref.  (1)). 
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It  is  clear  from  the  data  In  figures  1  to  7  and  fros  the 
data  in  table  4  that  transition  not  only  occurred  at  values 
of  Rea  less  than  Rea,c  also  occurred  where  Re^  was 
infinite.  Moreover,  except  for  the  29®  hemisphere-cone  at 
Mqo  -  2.32  and  at  2.47,  transition,  when  it  occurred, 

occurred  even  though  Reg  f,  was  much  larger  than  Reg  for  the 
entire  region  between  th4  stagnation  point  and  the  transition 
region.  Consequently,  it  is  concluded  that  for  laminar  flow 
to  exist  it  is  not  sufficient  that  Reg  be  less  than  Reg  c* 
Moreover,  it  is  also  concluded  that  the  probable  reason 'for 
the  occurrence  of  transition  even  though  Reg  was  less  than 
Reg  Cl  that  transition  was  not  caused  by  the  growth  of 
the 'small  two-dimensional  wave-like  disturbances  imposed  at 
one  Instant  as  assumed  in  the  theory  for  the  calculation  of 
Reg  c*  Because  the  values  of  Reg  are  approximate,  they  may 
be  ioo  large.  It  is  believed,  however,  that  they  are  not 
sufficiently  greater  than  the  correct  values  to  affect  these 
conclusions. 

Comparison  with  Calculated  Minimum  Transition  Reynolds  Numbers 


When  the  disturbances  to  the  laminar  flow  are  large  enough, 
transition  can  occur  even  though  Reg  is  less  than  Reg  c  (sss 
for  example,  ref.  (2)  and  ref.  (23)).  An  attempt  was' previously 
made  to  estimate  the  smallest  Reynolds  number  at  which  transition 
can  begin  (ref.  (2))  and  led  to  the  result  that  this  Reynolds 
number,  called  Reg  bi>  1^  the  value  of  Reg  at  which  the  local 
laminar  and  turbulent  friction  coefficients  are  equal.  Conse¬ 
quently,  in  order  to  calculate  the  value  of  Reg^m  at  an 
arbitrary  value  of  Mg  and  it  is  first  necessary  to  calcu¬ 
late  the  effect  of  Mg  and  on  the  laminar  and  turbulent 

friction  coefficients.  This  effect  was  calculated  by  use  of 
the  reference  enthalpy  method  (ref.  (24)).  Because,  however, 
of  some  uncertainty  in  the  ability  of  the  reference  enthalpy 
method  to  predict  the  turbulent  friction  coefficient  for  cold 
walls.  Regain  was  calculated  under  the  assumption  that  a  good 
estimate  for  Reg  q  is  the  value  for  an  insulated  wall,  even 

if  the  wall  is  colder  (see  ref.  (2)).  The  use  of  the  refer¬ 

ence  enthalpy  method  causes  a  decrease  in  the  estimated  value 
of  Reg  10  as  the  wall  temperature  is  lowered.  For  example, 
the  value  of  Reg  „  at  x  -  .7204  on  the  body  designated  as 
"spherical-segment-nose"  cylinder,  where  t^/To  is  about  .27 
and  Mg  is  4.064,  is  slightly  less  than  50,  whereas  the  value 
given  in  table  4  for  an  insulated  wall  and  shown  in  figure  6a 
is  95.4.  Consequently,  for  the  bodies  with  cold  walls  the 
values  of  Reg  „  given  in  table  4  may  be  as  much  as  double 
the  values  th&t  would  have  been  calculated  for  the  actual  wall 
temperature  ratios.  This  fact,  however,  does  not  change  the 
conclusion  Illustrated  in  figure  8,  namely,  that  the  values 
of  Reg  at  transition,  called  Reg^T>  were  in  all  cases  much 
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larger  than  the  values  of  Re^  ^  at  transition.  The  value  of 
Ree  'f  Is  taken  at  the  last  thermocouple  at  which  the  flow  was 
laminar.  Therefore,  the  transition  point  values  of  Re^  are 
actually  somewhat  larger  than  those  given  in  figure  8.  Because 
all  the  values  of  Red^'P  were  larger  than  the  values  of  Ree  ^ 
at  transition,  the  present  data  show  no  contradiction  of  t6e 
concept  or  of  the  method  for  estimating  Re^^n. 

Roughness  Reynolds  Number 

The  roughness  Reynolds  number,  Re^^,  was  calculated  because 
it  is  known  to  be  a  significant  measure  of  the  ability  of 
roughness  to  cause  transition.  For  example,  investigators 
have  found  that  the  average  value  of  Rej^  needed  to  produce 
transition  at  a  1/4-inch  wide  strip  of  sandpaper-type  roughness 
in  subsonic  flew  is  about  400  (see  for  example  ref.  (25)). 
Moreover,  reference  (25)  states  that  the  transition  Reynolds 
number  based  on  x,  decreases  to  95  percent  of  its  value  without 
roughness  when  Re^^  for  a  1/4-lnch  wide  strip  of  sandpaper  is 
between  178  and  330. 

For  the  bodies  of  revolution  considered  in  the  present 
investigation,  Re^^  begins  at  zero  at  the  stagnation  point 
and  then  Increases;  it  may  have  one  or  more  local  maxima  on 
the  body.  For  the  cases  analyzed  in  the  present  investigation 
the  maximum  values  of  Re|^  fall  between  .0082,  calculated  for 
the  29<>  hemisphere-cone  at  2.32  and  the  value  30.80 

calculated  for  section  ALM  of  the  Elliptical-Nose  Cylinder. 

These  values  of  Reic  are  much  smaller  than  the  values  of  RSk 
for  which  strips  of  sandpaper -type  roughness  first  affect 
transition.  For  example,  the  largest  value,  30.80,  is  about 
l/6th  of  the  smallest  value,  178,  reported  for  a  1/4-inch 
wide  strip  of  sandpaper  (ref.  (25)). 

In  order  to  see  whether  or  not  there  is  any  relation 
between  the  largest  value  of  Re^  in  the  laminar  boundary  layer 
on  a  body  of  revolution  and  Ren  j,  the  values  of  Re^^T  ssre 
plotted  against  the  largest  value  of  Re^  ahead  of  transition 
for  five  of  the  seven  cases  analyzed  (fig.  9).  (See  also  Table  5) 
Two  of  the  cases  do  not  appear  on  figure  9  because  for  one  of 
them,  the  Hemisphere-Cylinder,  the  flow  was  already  turbulent 
at  the  first  thermocouple  located  at  11  1/2^  from  the  nose.  Con¬ 
sequently,  the  last  and  only  thermocouple  that  indicated  laminar 
flow  was  the  one  at  the  nose.  The  original  investigators 
(ref.  (11))  suggest  that  perhaps  irregularities  other  than  the 
five-microinch  surface  roughness  caused  the  early  transition. 

The  other  case  of  the  seven  that  does  not  appear  in  figure  9 
is  that  for  the  Spherical-Segment-Nose  Cylinder.  In  this  case 
transition  did  not  occur  on  the  portion  of  the  body  containing 
thermocouples  and  so  values  of  Re^  j  could  not  be  calculated. 
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The  maximum  value  of  Rej^  ahead  of  transition,  called  Re^  K, 
rather  than  its  value  at  transition,  was  used  because  for  a  ’ 
surface  with  continuously  distributed  roughness,  Rei^  H  rather 
than  Rej^  j,  is  a  measure  of  the  largest  disturbance  introduced 
into  the 'boundary  layer  by  the  roughness  and  so  is  probably  a 
better  measure  of  the  effect  of  roughness  on  transition  than 
Rej^  T*  analogous  to  the  use  of  Rej^  when 

transition  is  caused  by  a  roughness  strip.  Thus,  when  tran¬ 
sition  occurs  downstream  of  the  roughness  strip,  Rei^  t  zero 
and  the  relation  is  between  Ree^T  and  Rejc  M«  Actually,  a 
better  criterion  than  either  Rej^  n  or  ReK,T  is  probably  one 
that  depends  on  the  distribution 'of  Re^^  in  the  entire  region 
upstream  of  transition. 

It  may  be  of  Interest  to  note  that  the  data  for  a  surface 
completely  covered  with  sand  roughness  (ref.  (26)  and  page 
450  of  ref.  (22))  show  that_the  transition  point  is  first 
affected  by  roughness  when  reaches  about  120,  or  about 

60  when  the  nominal  height  Ts.  is  substituted  for  the  equivalent 
sand  roughness  Eg.  The  associated  value  of  Rex  T  is  about 
.664  X  10°.  For  these  data  transition  is  therefore  first 
affected  when  Rek,ll  is  near  60.  On  the  other  band,  a  calcu¬ 
lation  of  Re^  7  for  these  data  results  in  a  value  of  about 
1.5.  This  number  is  obtained  from  the  definition 


V 


where 


and  the  relation. 


(incompressible  flow) 
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Then 


but 


[jj  J  ^ 


for  the  zero  pressure  gradient  data  of  reference  (26) . 
Therefore 
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For 


/  ACx  \  y  J 


3nd 


=■  Co 


the  result  is 

1-^ 

That  Is,  the  value  of  Rek>T  when  transition  first  begins  to 
move  forward  on  a  surface  with  continously  distributed  roughness 
is  much  less  than  values  of  Re^  previously  quoted  for  the  first 
effect  of  roughness  strips  on  transition  (ref.  (25)).  Note, 
however,  that  values  of  Rek,T  not  far  from  1.47  have  previously 
appeared  in  the  literature.  For  example,  in  reference  (27)  a 
value  of  Re^  <p  as  low  as  7.35  is  given.  Moreover,  In  reference 
(19)  a  still’ lower  value  of  Re^  namely,  about  two  is  given. 

In  this  case  the  value  of  Re^  ai  the  location  where  7  is  equal 
to  7  is  about  400;  this  value  is  thus  equal  to  Re^  ii>  Conse¬ 
quently,  it  appears  that  for  surfaces  completely  covered  with 
roughness  the  value  of  Re]^  X  can  be  near  unity  where  transition 
is  first  affected  by  the  roughness.  The  smallest  value  of 

previously  given  in  the  literature,  seems,  however,  to 
be  that  of  reference  (26),  that  is,  a  value  near  60. 


The  data  in  figure  9  cover  a  range  of  Re^  m  ^rom  about 
.003  to  about  31,  that  is,  from  about  1/7000  io  about  1/2  of 
the  value  60.  The  largest  value,  30.8,  occurred  on  a  patch 
of  roughness  that  caused  Reo  x  drop  from  about  1082  to  about 
177  on  the  l%.lliptical-Nose  Cylinder.  Consequently,  transition 
is  affected  when  Re^  g  is  as  low  as  30.8.  Moreover,  the  data 
in  figure  9  indicate ’ that  the  two  points  for  the  Elliptical- 
Nose  Cylinder,  namely,  the  points  for  ReK^y  of  30.8  and  .01825, 
are  consistent  with  the  variation  between  Reg  x  M 

indicated  by  the  other  four  sets  of  data.  It’ is  emphasized 
that  all  the  data  of  figure  9  are  for  blunt  bodies  of  revolution 
in  supersonic  flight  with  Reg^^.  much  larger  than  Reg. 

The  line  drawn  through  the  data  in  figure  9  is  a  least- 
squares  line  for  log  Reg^X  SEkinst  log  Reg^lf.  Its  equation  is 

./8» 
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Although  the  equation  of  the  line  Is  given,  it  is  not  concluded 
that  figure  9  establishes  a  connection  between  Reg  7  and  Rej^  l|. 
The  reason  is  that  an  application  of  the  Information  in  refer¬ 
ence  (26)  shows  that  the  scatter  of  the  data  around  the  least- 
squares  line  is  too  large  to  allow  the  conclusion  to  be  drawn 
from  only  five  sets  of  data  that  figure  9  establishes  a  re¬ 
lation  between  Reg  x  and  Re)(  n  for  blunt  bodies  of  revolution 
in  supersonic  flight  with  Reg,c  much  greater  than  Reg. 

On  the  other  hand,  even  if  there  really  were  a  relation 
between  Reg  j  and  Re^^  y  somewhat  like  that  shown  by  the  least- 
squares  Hue  of  figure  9,  some  scatter  of  the  data  would  still 
be  expected.  The  reasons  are:  first,  the  values  of  Reg^T 
were  computed  for  the  rearmost  thermocouple  at  which  the  flow 
was  laminar.  Because  transition  often  began  between  two 
thermocouples,  its  precise  location  and,  consequently,  the 
precise  value  of  Reg  j,  could  not  be  found.  Second,  the 
value  of  the  surf ace ’ roughness  was  not  obtained  in  the  same 
way  in  all  the  tests  (see  Table  3).  Consequently,  some  of 
the  roughness  Reynolds  numbers  probably  are  not  directly 
comparable.  Note  that  Re^^  is  proportional  to  the  roughness 
height  squared.  Third,  the  surfaces  were  not  all  finished 
by  the  same  process  and  so  two  surfaces  with  the  same  measured 
roughness  height  can  affect  transition  differently  (see  ref. 
(29)).  The  fourth  reason  is  that  not  every  part  of  the  surface 
of  the  bodies  could  be  examined  with  a  microscope.  Fifth, 
these  roughness  data  are  for  the  bodies  before  flight.  Because 
the  bodies  could  not  be  examined  after  flight  the  conditions 
of  their  surfaces  when  transition  occurred  are  not  really  known. 
These  reasons  would  probably  be  sufficient  to  cause  a  fairly 
large  amount  of  scatter  even  if  there  were  a  strong  connection 
between  Reg  x  and  Re^^  y.  In  spite  of  these  reasons  for  the 
data  to  scatter,  it  is  probably  not  correct  to  conclude 
anything  more  from  figure  9  than  that  the  previous  impression, 
namely,  "the  smoother  the  surface  the  greater  the  likelihood 
of  extensive  regions  of  laminar  flow,"  is  re-enforced  by  these 
data. 


It  is  remarked  that  a  plot  of  Reg  x  against  (^w^e^T>  shown 
in  figure  10,  (see  also  Table  5)  Indicates  a  very  much  weaker 
connection  between  Reg  x  snd  (t^/tc)T  than  Indicated  in  figure  9 
between  Reg  x  K^k  M*  Note  the  large  amount  of  scatter.  The 
line  drawn  in_figure  lO  is  a  least-squares  line  for  log  Reg  ^ 
against  log  (tw/te)x  snd  has  the  equation  ’ 

=  ^5-0 
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SUMMARY 

Because  of  the  Importance  of  aerodynamic  heating  at  high 
speeds  and  because  the  rate  of  heat  transfer  is  much  greater 
for  a  turbulent  than  for  a  laminar  boundary  layer  it  is 
important  to  be  able  to  estimate  the  location  at  which  the 
flow  changes  from  laminar  to  turbulent.  The  main  result  of 
previous  theoretical  work  is  a  theory  that  can  be  used  to 
estimate  the  location  at  which  the  laminar  boundary  layer 
becomes  unstable  once  the  necessary  parameters  of  the  laminar 
boundary  layer  are  known.  In  the  present  investigation  these 
parameters  are  calculated  by  the  Cohen  and  Reshotko  method  for 
seven  blunt  bodies  of  revolution  flying  at  supersonic  speeds. 

A  comparison  of  the  results  of  the  stability  theory  with 
transition  points  determined  from  heat-transfer  distributions 
obtained  for  these  bodies  by  previous  investigators  then  shows 
that  when  transition  occurred  it  took  place  even  though  the 
boundary  layer  was  computed  to  be  very  stable  for  the  entire 
region  between  the  stagnation  and  the  transition  point.  It 
is  suggested  that  in  these  cases  transition  probably  was  not 
caused  by  the  growth  of  the  small  two-dimensional  wave-like 
disturbances  imposed  at  one  instant  as  assumed  in  the  stability 
theory  used  in  the  present  investigation. 

In  every  case  transition  occurred  at  a  larger  boundary- 
layer  Reynolds  number  than  the  estimated  minimum  Reynolds 
number  for  transition.  Consequently,  no  contradiction  of 
the  assumption  that  there  is  a  minimum  transition  Reynolds 
number  and  no  disagreement  with  the  results  of  the  method  for 
estimating  this  Reynolds  number  was  found. 

A  plot  of  the  calculated  boundary-layer  Reynolds  number 
at  transition  against  the  calculated  maximum  roughness  Reynolds 
number  ahead  of  transition  for  the  five  of  the  seven  sets  of 
data  for  which  transition  data  were  available  apparently 
showed  a  connection.  A  further  examination  indicated,  however, 
that  the  scatter  of  these  data  is  too  large  to  allow  the 
conclusion  to  be  drawn  from  only  five  sets  of  data  that  there 
really  is  a  connection.  The  conclusion  seems  to  be  only  that 
the  likelihood  of  obtaining  large  regions  of  laminar  flow 
increases  as  the  body  is  made  smoother.  It  is  remarked  that 
the  dependence  of  the  boundary-layer  Reynolds  number  at 
transition  on  the  wall  temperature  ratio  at  transition  was 
very  much  weaker  than  on  the  maximum  roughness  Reynolds 
number • 
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BODY  PRESSURE  DlStR IBUTION ,  WALL  TEMPERATURE  DISTRIBUTION. 
AND  RADIUS  DISTRIBUTION  DATA 


Bod)  Pressure  Distribution 

Wall  Temperature  Distribution  Radius  Distribution 
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Table  5 


CALCULATED  VALUES  OF  Re^  j,.  AND  {  — 


FOR 


le/T 


FIVE  OF  THE  SEVEN  BLUNT  BODIES  OF  REVOLUTION 


®®e,T 

IW\ 

Vte/  T 

1557 

.008195 

1.011 

1555 

.009850 

.9634 

888 

.01115 

.7802 

713 

.01243 

.6776 

632 

.01426 

.6282 

642 

.01579 

.5959 

251 

4.270 

.6934 

274 

6.600 

.6038 

293 

9.060 

.5472 

297 

10.06 

.5108 

303 

11.02 

.4838 

579 

.9400 

.8148 

449 

.9840 

.7143 

465 

1.144 

.6516 

519 

2.100 

.4988 

546 

2.717 

.4495 

278 

1.898 

.3542 

277 

1.790 

.3949 

304 

1.055 

.1675 

397 

1.610 

.1464 

446 

1.710 

.1417 

1082 

.01825 

.2271 

177 

30.80 

.1866 
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(b)  LOCAL  MACH  NUMBER,  LOCAL  WALL  TEMPERATURE  RATIO 
AND  LOCAL  ROUGHNESS  REYNOLDS  NUMBER 


THERMOCOUPLE  LOCATIONS 
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(a)  BOUNDARY  LAYER  REYNOLDS  NUMBERS 

FIG.  3  HEMISPHERE-CYLINDER  IN  FLIGHT  AT  A  MACH  NUMBER  OF 
5.50  AND  A  REYNOLDS  NUMBER,  Re  ,  OF  9.75  x  10® 
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(b)  LOCAL  MACH  NUMBER,  LOCAL  WALL  TEMPERATURE  RATIO,  AND 
LOCAL  ROUGHNESS  REYNOLDS  NUMBER 
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(b)  LOCAL  MACH  NUMBER,  LOCAL  WALL  TEMPERATURE  RATIO, 
AND  LOCAL  REYNOLDS  NUMBER 
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(b)  LOCAL  MAC.  NUMBER,  LOCAL  WALL  TEMPERATURE  RATIO. 
AND  LOCAL  ROUGHNESS  REYNOLDS  NUMBER 
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FIG.  9  TRANSITION  REYNOLDS  NUMBER  RCo  AND  MAXIMUM  ROUGHNESS 
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